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PREFACE 


Objective 

The purpose of this study is to achieve a preliminary 
design for a shuttle based BUV imaging spectrometer 
that meets the mission requirements. 

Scope of Work 

Using the present BUV instrument aboard NIMBUS 4 
as a baseline, compare point scanner mechanisms 
and spatial multiplex scanning systems based on 
sensitivity, field of view and simplicity. Include in 
the comparison spatial and spectral scanning and 
multiplexing techniques, with a view at optimizing 
the overall system for the basic performance 
requirements. Select an optimum sensor configuration 
and perform a conceptual design study for the selected 
sensor. 

Conclusions 

An optimum sensor configuration using a "pushbroom" 
spatial scanner with a Hadamard spectral scan was 
selected. This system meets the requirements for a 
shuttle based mission. A conceptual design was 
performed on the selected sensor and monochromator 
ray traces, mechanical layouts and the electronic 
block diagram were generated. 

Summary of Recommendations 

A High Resolution Ozone Mapper should be bread- 
boarded and tested. 
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INTRODUCTION 


1 . 0 

This final report on the High Resolution Ozone Mapper (HROM) 
is submitted as required under Item 9 of the Statement of Work 
for Contract NAS5-22326. 

The purpose of the study was to achieve a preliminary design 
for a shuttle based HROM Imaging Spectrometer System that 
meets the requirements of the subject work statement. This 
report details the accomplishments of the study. 

The performance requirements of the work statement constitute 
a significant increase in sensitivity over the Backscatter 
Ultraviolet (BUV) Instrument aboard NIMBUS 4. This instrument 
was used as a baseline in the study. The BUV instrument for its 
size and configuration represents a limit in performance for 
state-of-the art technology. Thus, in order to realize the 
objective of the study, a different type of approach 
was required. A more satisfactory approach was 
realized which combined the obvious advantages of a "pushbroom" 
scanner using a linear array of detectors for spatial scanning 
with the improvements derived from using a multiplex system for 
spectral scanning. This report contains the rationale, analysis 
and preliminary mechanization design for an instrument which 
uses this approach and meets the objective of the work state- 
ment. 

Section 2 of the report discusses the requirements imposed 
on the instrument to insure faithful reproduction of the ozone 
profile and radiance level as well as a listing of the overall 
system requirements. 
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Section 3 of the report details the results of the candidate 
sensor tradeoff. In this section, the types of systems initially 
considered are enumerated. Then, the rationale and analysis 
by which the number of systems were reduced is detailed. 

Finally, a candidate tradeoff matrix is generated and the final 
candidate systems (2) selected. 

Section 4 contains the performance evaluation study results. 

The section briefly discusses the fact that the conventional 
system can be evaluated by extrapolation from the BUV instrument, 
but the Hadamard System must be discussed in detail. Thus, 
the section mainly discusses Hadamard Spectrometers. The 
sources of noise and the effects of these noise sources on the 
sensor are delineated. Also, the section discusses errors and 
their correction. Next, a preliminary design of the Hadamard 
mask is generated. Finally a comparison of the two candidate 
sensors is presented, using signal-to-noise ratio as a criteria. 

Section 5 contains the selection of the recommended system. 

The section briefly recounts the rationale used to select the 
final candidate system. 

Section 6 contains the results of the comceptual design activity 
task. In this section, a preliminary system mechanization is 
presented. The section contains the mechanical, optical and 
electronic design for an instrument which meets the requirements 
for the mission. 

This report has been prepared by AS&E, incorporating contributions 
in the Requirements, Mechanical Design and Optical Design 
Section from Beckman Instrument, Inc. , and in the Hadamard 
Multiplexing Section, from Spectral Imaging, Inc. 
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2. 0 EXPERIMENT REQUIREMENTS 

The High Resolution Ozone Mapper (HROM) System will provide 
the ultraviolet sounding and mapping data necessary to calculate 
and map the changing distribution of the ozone layer surrounding 
the earth. 

To perform these functions, analysis shows that very high spectral 
resolution must be attained and the stray light reduction and 
electronic dynamic range must be large. 

2. 1 Spectral Requirements 

o 

The spectral bandpass must be narrow (~ lOA), the wavelength 

accuracy high (-0, 2 to 0. 5^) and the dynamic range very large 
5 7 

(10 to 10 ) because of the steepness of the ozone 
absorptive edge. The successful performance of the Nimbus 4 
BUV for four years in space shows that these requirements 
can be met reliably. 

The extremely rapid wavelength dependence of the earth's 
ultraviolet albedo at the edge of the strongly absorbing 
ozone band determines the instrument requirements for 
spectral bandpass, wavelength repeatability, dynamic 
range, and stray light. 

One can see why the instrument spectral selection requirements 
need to be as tight as they are by considering the wavelength 
variation of the received albedo radiance. A plot of received 
radiance as a function of wavelength is shown in Figure 2-1. 

This plot is taken from Nimbus 4 BUV data reported by 
Heath and Heaney. The radiance which would be observed if the 
solar energy were diffusely reflected from a perfect diffuser 
with 100% albedo is shown for comparison. The ozone 
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2500 3000 35 
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Fig. 2-1 Observed Radiance Versus Wavelength (BUV Orbit 100) 




absorption causes a decrease in observed radiance at a 
rate of 2 decades per 120 % in the neighborhood of 3050 % 

To infer the ozone profile reliably from the data, the atmospheric 
albedo must be determined within limits of 1 to 3 percent, 
which requires measurements of the radiance to 1 percent 
or better. The bandpass, wavelength accuracy and 
repeatability of the spectrometer must be so chosen as to 
contribute only a small part of this error. 

2.2 Wavelengtn Accuracy and Repeatability 

In the neighborhood of 3050 % the radiance (N) changes at a rate 

o 

of two decades over 120 A (or a factor of e = 2. 718 every 
26 2o) The fractional error in radiance resulting from a 
wavelength measurement error A x is then: 

AN = AX (2. 1) 

N Xo 

where Xq “ 26%, To obtain an error of ll percent, the value 
of A X must be less than 0. 26 %. 

A slightly reduced value for absolute accuracy of the wavelength 
can be accepted because several measurements are made but the 
repeatability of the wavelength to -0. 2 to 0. 5^ is a necessary 
requirement, 

2. 3 Bandpass 

If the instrument bandpass is too wide, the profile will be 
smoothed with loss of information. For a spectrum varying 
over a limited range, x to x the radiance N is given by the 
expression: 
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(x- Xl)/Xo 


( 2 . 2 ) 


N = N (Xj) e 
where N (X^) = the radiance at X^ 

The average radiance over a triangular bandpass of width at 
half-height aX centered at X^ is: 


N = 




(2.3) 


This average value deviates from N { Xj^) by an error which, 
to second order in A xAq given by 

_ 2 

-,4 N . = — . (2.4) 

N(Xi) 12\^Xo/ 

o o 

For aX = 10 A, Xq = 26 A, the shift in the radiance is 1. 2 
percent. Although the effect can be corrected by deconvolution, 
it is desirable to keep it as small as possible. This consideration 

explains why the spectral bandpass must be as small as 10 A. 

2, 4 Resolution 

The inherent resolution of the instrument must be better than 
o 

1 A (Rayleigh criterion). If the inherent resolution, which 
depends on wavelength, were larger, the overall bandpass would 
exhibit wavelength-dependent variations larger then 1 percent, 
leading to a comparable error in the spectral radiance. 

2.5 Stray Light 

Stray light refers to the light of the wrong color detected by the 
spectrophotometer. (For a filter radiometer, this is conven- 
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tionally called "out of band response, ") Only a double grating 

0 

monochromator can provide the stray light rejection of 10 
required to meet the dynamic range requirements. 

2. 6 System Specification 

Table 2. 1 contains a list of the system preliminary specifications 
as contained in the Statement of Work. 
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TABLE 2o 1 


PRELIMINARY SPECIFICATIONS 


Orbital Altitude 

400 Km 

Total Field of View 

o o 

7 X 104 

o O 

Instantaneous Field of View 

7x7 

Number of Detectors 

15 

Wavelength of Spectral Bands 

2555, 2735, 2830, 
2876, 2922, 2975, 
3019, 3058, 3125, 
3175, 3312, 
3398A® (12 Bands) 

Band pass (Half Power) 

o 

lOA 

Spectral Resolution 

o 

1 A 

Stray Light Rejection 

10® 

Polarization Sensitivity 

<5% 

Signal Dynamic Range 

3 

10 

-4 

S/N Ratio 

30 @ 2.^5 X 10 Q 
Erg/cm'^ Ster. A g 

Radiometric Precision 

1% 

Overlap / Underlap 

<5% 
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CANDIDATE SYSTEM TRADEOFF ANALYSIS 


3. 0 

3, 1 Types of Systems Considered 

Using the assumption that each detector in a multielement 

system must view all the spectral regions (i. e., separate 

detectors cannot be used for each spectral region), the 

candidate system contained in Table 3-1 were generated. A 

description of each of the systems follows: 

3. 1, 1. Candidate Sensors Which View the Entire 104 ° 

o 

by 104 Field of View 

a) Hadamard Spatial and Hadamard Spectral Multi- 
plexing Sensor 

Light from the total field of view is collected by 
the aperture and focused onto a two dimensional 
mask with square open and blocking elements 
(looks like a cross-word puzzle). This mask 
rotates (or steps) at a rate commensurate with the 
frame time (for a 109 sec, frame time each mask 
position must be held for 35 msec). In each 
mask position 255 (15 x 17 elements) of the 
elements are displayed through an area determined 
by a rectangular field mask. One hundred twenty 
eight (128) of these elements pass light, and one 
hundred twenty seven (127) elements block it in 
any given mask position. There are two hundred 
fifty-five (255) Independent positions for this mask. 
The mask contains a total of 957 elements for 
each frame scanned. 
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TABLE 3-1 


TYPE OF SYSTEMS CONSIDERED 


1 . 


2 . 


View Entire 104° FOV 



View FOV of 104° by 7X° 



3. View Hor. Line 104° Wide by 
7° High 


/|\ 


4. View Vert. Line 104 High by 
7° High 


5. View One Spatial Res Elm 
70 by 70 


Where 


Spatial 

a. 2 Dim H 

b. Image Plane Scanner 

c. Image Plane Scanner 

d . Image Tube 

e. Mosaic Array 


a . 1 Dim H 

b. Image Plane Scanner 

c. Image Plane Scanner 

d. Image Tube 

e. Mosaic Array 


a. 1 Dim H 

b. 1 Dim H 

c. 15 element lin array 

d. 15 element lin array 


a . 1 Dim H + OPS 

b. 1 Dim H + OPS 

c . Linear Array + OPS 

d . Linear Array + OPS 


a. OPS 

b. OPS 


2 Dim H Two dimension Hadamard Mask 


1 Dim H One dimension Hadamard Mask 


C Same Spectral Scanning as BUV 

OPS Object Plane Scanner 


Spectral 

1 Dim H 
1 Dim H 
C 

Undecided 

Undecided 


1 Dim H 
1 Dim H 
C 

Undecided 

Undecided 


1 Dim H 
C 

1 Dim H 
C 


1 Dim H 
C 

1 Dim H 
C 


1 Dim H 
C 
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After passage through the mask the modulated radiation 
is dispersed by the monochromator and refocused onto 
a second mask. This second modulating mask is a one 
dimensional spectral mask. It consists of 23 slits, some 
open, some opaque. Twelve of these slits are exposed 
at one time. There are also twelve independent positions 
for this mask. Each position must be held for a full cycle 
of the first mask (approximately 9 sec.). After passing 
through the second mask the spatial and spectral encoded 
radiation is de-dispersed and focused onto a PMT. 

The electrical output of the PMT is conditioned and con- 
verted to digital format. The digital words are then 
telemetered to the ground station, where the Hadamard 
transform is performed on the data to reconstruct the 
scene. 

b) Image Plane Scanner Spatial Scan With A Hadamard 
Spectral Scan 

An image of the entire 104° square field is produced using 
a lens system and an image plane scanner is used to 
interrogate this field pixel by pixel. A modified BUV 
monochromator is employed with a single phototube and 
a Hadamard one dimensional spectral mask scans the 
wavelengths . 

Several image plane scanner types could be used, but in 
a system with an aperture less than approximately 8 inches 
there is no advantage in mechanical complexity over an 
object plane scanner. The simpler types of image plane 
scanners counter rotating wedges, reticles would require 
a form of image motion compensation to direct the light 
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through the monochromator. The more complex image 
plane scanners such as the squirrel cage or the roof 
mirror are significantly more complex than an object 
plane scanner such as that used in TOMS and would be 
difficult to build in the small size required by HROM, 
External reflecting polygons and two axis oscillating 
mirrors are about equal in complexity to the TOMS 
scanner and would have the same system throughput. 

The radiation from each pixel in the image plane is 
dispersed by the monochromator and focused onto a 
Hadamard one dimensional mask containing 23 slits. 

Twelve slits are exposed at one time and there are 
twelve independent positions. Each position must be 
held for one complete scan of the image plane (one 
field) approximately 9 sec. After passing through the 
mask the modulated radiation is de-dispersed and focused 
onto a PMT. 

The electrical output from the PMT is conditioned digitized, 
telemetered and Hadamard transformed to reconstruct the 
original scene . 

c) Image Plane Scanner Spatial Scan With Conventional 
Monochromator 

This system is the same as (b) above except for spectral 
scanning. The present BUV monochromator is employed with 
a single phototube and a step scan through the wavelengths. 
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d) Image Tube Systems 

Ultraviolet Image tubes are available with sensitivity 
adequate to meet this specification. One such tube is 
the Digicon, wherein photoelectrons from an end-on 
semi-transparent photocathode are accelerated and mag- 
netically focused on an array of silicon diodes. Tubes 
of this type with 1024 elements have been built. 

While it is not clear how the spectral requirements 
could be met, the extreme sensitivity of these devices 
makes the approach interesting. 

A system built recently at the University of Texas was 
capable of counting as few as one photoelectron every 
four minutes. Unfortunately, saturation of these devices 
occurs at a total count of about 6000 photoelectrons. To 
accommodate the dynamic range of HROM the system in- 
tegration time would have to change with signal level by 
about three orders of magnitude. While circuitry could be 
designed for this, the low saturation level corresponds to 
a limitation on measurement accuracy of about 2%. Another 
source of error arises from the fact that the diode array is 
also sensitive to light from 4000 to 10,000 S. In the 
University of Texas system, which was operating at approx- 
imately 4000 S, the error from a Tungston source was 2%. 

For HROM it would probably be much larger. Two additional 
problems, radiometric nonlinearity and the necessity for 
refrigeration to approximately -80°C reduce the attractive- 
ness of this system in spite of its sensitivity. 
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e) 


Mosaic Arrays 


Mosaic arrays are available with sensitivity adequate 
to meet this specification. The arrays are presently 
silicon detectors with CCD or CID readouts. Arrays 
of this type have been produced with up to 1024 
elements. 

The major drawback to this type of system is in the 
dynamic range. Present devices only operate over 
two to three orders of magnitude. The major problem 
area is saturation in the CCD device. Thus while 
mosaic arrays are attractive from a system mechaniza- 
tion standpoint, the devices are not sufficiently de- 
veloped for this application. 

3,1,2. Candidate Sensors which View A Portion of the Total 
Field of View 

a) Hadamard Spatial and Hadamard Spectral Multi- 
plexing Sensor 

This sensor is similar to sensor la, the only difference 
is the two dimensional spatial mask, which in this 
case is 15 elements wide by x elements high (depend- 
ding on the vertical field of view). 

b) Image Plane Scanner For Spatial Scan With a Hadamard 
Spectral Scan 

It is possible to construct an imaging system which 
would view a total field of 104° in one dimension 
and less than 104° in the other dimension. Both 
systems Ic) and Id) above could be configured in 
this manner but there is no obvious advantage other 
than minimizing the effects of sun angle changing 
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during the period of measurement. The system 
mechanization considerations are the same. 

c) Image Plane Scanner For Spatial Scan With A 
Conventional Monochrometer 

See Section 3. l.l.c. 

d) Image Tubes 

This system is similar to system Id. It is possible 
to construct an imaging system which only views 
a portion of the field. However, there is no obvious 
advantage to this system. 

e) Mosaic Array 

This sensor is similar to sensor ie. , the only difference 

is in the number of detectors in the array. 

3. 1, 3, Pushbroom System (View a line 104° wide by 7° high) 

a) Hadamard Spatial and Hadamard Spectral Multiplexing 

Sensor 

The light is collected by the aperture slit and dis- 
persed perpendicular to the 104° direction by the 
monochromator and focused onto a two dimens ioxial 
Hadamard mask. The mask encodes the spatial informa- 
tion in the direction parallel to the 104° and spec- 
trally in the direction of dispersion. The mask will expose 
15 X 12 elements for each position and have 15 positions. 
Possibly the mask will be on a wheel with 12 one dimen- 
sional 15 element separate masks which rotate into 
place. The energy passing through the mask will be 
de-dispersed and focused onto a single PMT. 
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b) Hadamard Spatial and Conventional Spectral Sensor 

The light is collected by the aperture and focused 
onto a one dimensional Hadamard mask containing 
29 slits. The mask exposes 15 slits at any one time 
and has 15 independent positions. The encoded 
radiation passing through the mask is dispersed by a 
monochromator and focused on an exit slit. The mono- 
chromator is a modified BUV design with a single 
phototube and a step scan through the wavelengths. 

c) Linear Array Spatial Scan with a Hadamard Spectral 
Scan 

A lens is used to image a field 104° wide and 7° high. 

A linear array of 15 photomultipliers detects the energy 
in the 15 pixel wide field. Since system sensitivity 
increases with the square root of the number of detec- 
tors, an improvement in signal-to-noise ratio of four- 
fold is achieved. The monochromator is a modified BUV 

* 

design where the dispersed radiation is focused onto 
a one dimensional Hadamard mask containing 29 slits. 
Twelve of the transmissive and opaque slits are ex- 
posed to the radiation for a time comenserate with 
the frame time. There are 12 independent mask positions. 
The encoded radiation passing through the mask is 
de-dispersed and focused onto the linear array of 15 PMT's. 
The linear array does the spatial scanning. 

d) Linear Array Spatial Scan with a Conventional Mono- 
chromator 

This system is the same as 3c except the present BUV 
monochromator concept is used to provide the spectral 
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information and the only modification is a rotation 
of the monochromator so that dispersion is in the 
direction of the flight vector. The system is basically 
the same as the present BUV except for the imaging 
lens and the increased number of detectors and pre- 
amplifiers. 

It would also be possible to configure the system so that 
the extra detectors are used in the spectral domain. The 
matching of the photomultipliers would then become 
a major instrument design problem. Gain stability, 
linearity over the total dynamic range and space radi- 
ation effects could lead to errors unable to be removed by 
calibration. For this reason no systems were considered 
wherein multiple detectors were used for the different 
ozone absorption wavelengths. 

3, 1.4. Object Plane Scanning Systems (Views a Reid of View 
104° high by 7° wide) 

An object plane scanner could be employed with an 
array of detectors. The system sensitivity would 
be improved by the square root of the number of detec- 
tors used. The system telescope requires an imaging 
lens in addition to the mechanical scanner. If 15 
detectors are used there would be no advantage to these 
systems over the pushbroom scanners described above 
and they would have the additional complexity of the 
added scanning mechanism. If less than 15 detectors 
are used, the sensitivity improvement would be diminished. 
One advantage of an object plane scanner over a push- 
broom is that they are somewhat easier to calibrate. The 
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calibration can be done during the dark period of the 
scan cycle. With a pushbroom, the calibration 
sequence results in loss of ground coverage for the 
duration of the sequence. 

3, 1.5. Point Scanning System (View a Field of View 7° by 7°) 

An object plane scanner could be employed with a 
single detector. The scan is perpendicular to the 
flight vector. The system telescope requires a narrow 
field of view lens (7°) in addition to the mechanical 
scanner. For system 5a a Hadamard mask, which is the 
same as that used in system 3c, is used for spectral 
scanning. In system 5b, a BUV monochromator with 
a stepping through the wavelengths is used. This 
type of system is very similar to the present BUV 
and the sensitivity does not meet the requirements. 
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First Reduction of Candidate Systems 


3. 2 

Referring to Table 3-1 several systems were eliminated by 
inspection. These included: 

1. Image Plane Scanning Systems (Systems lb,lc, 2b and 2c). 
Since the aperture diameter for the HROM system is small, 
an image plane scanner has no size advantage over an 
object plan scanner. Also the image plane scanner 
requires wide angle optics. Therefore, it was decided to 
eliminate these systems. 

2. System Ulilizing Image Tubes (Systems Id and 2d) 

Since image tubes do not have the capability to meet the 
dynamic range requirements, these systems were also 
eliminated from further consideration. 

3. Systems Employing Mosaic Arrays (System le and 2e) 

Mosaic arrays of detectors with CCD readout were 

considered. However, state-of-the-art devices only 
4 

have a 10 dynamic range capability. Thus, these systems 
were eliminated. 

In this manner, the number of candidate sensors to be analyzed was 
reduced to 12 systems. 

3. 3 Signal-To-Noise Ratio Evaluation 

In order to further reduce the number of candidate systems, the 

systems were compared on a basis of signal-to-noise ratio. 

Equations were derived for the signal-to-noise ratio for both 
conventional and Hadamard scanning systems and these 
equations were used to calculate the SNR for each of the twelve 
systems. 
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3. 3. 1 


Conventional System ( point scanning system) 


For a photomultiplier system: 


S N 

— = s 



N + N^l 

1/2 

(3.1) 


s D 1 




where: 

N = Number of photoelectrons per second from the source 
s 

Af = System bandwidth 

Nj^ = The sum of the number of photoelectrons per second 
from the detector dark current, the trapped electron, 
X-rays and protons. 


3. 3. 2 Hadamard Multiplex System 

For a Hadamard sensor, since more than one spatial element is 
viewed simultaneously, equation (3. 1) becomes: 



where: 

M = Number of spatial elements scanned (Pixels) 

P = Number of spectral elements 

N = Number of masks required (transmission factor) 
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however: 



IS only the ratio for 1 field (one look) 



But in a full frame each pixel is viewed (M x N) times and 
since the S/N ratio improves by the square root of the number of 
looks 


/s\ = (s) ^ (M X P) (3. 4) 

Wh" Wh’ 

substituting eq. (3. 3) into eq. (3. 4): 



(M X P) 
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simplifying eg. (3. 6) 


(M K P) 


1/2 


H" 


2NAf 


1/2 


N /M X P\ . _ 


1/2 


(3. 7) 


For noise in signal limited operation equation (3. 7) becomes: 

1/2 


N 




(2NAf) 


1/2 


(Ns» Nd) 


(3.8) 


For additive noise limited operation eg. (3. 7) becomes: 

1/2 


(n)h„ 


N (M X P) 
s 


2N 


(Nd» ^s) 


(3.9) 


3. 3. 3 Additive Noise ( N^^) 

Figure 3-1 shows the S/N ratio for the twelve candidate 
systems. For this analysis the value of was set 

3 

at 2 X 10 photo electrons per second for a minimum value and 

7 

1. 95 X 10 photo electrons per second for the maximum value. 

The curves were generated using equations (3. 1) and (3. 7) 
where the values for the parameters are contained in Table 3. 2. 
The minimum value of additive noise (N^^) was obtained from the 
"quiet" values obtained from the BUV flight. The maximum value 
of additive was obtained from calculation at 400 Km in the 
South Atlantic Anomaly. 
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SIGNAL-TO-NOISE RATIO 
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SYSTEM IDENTIFICATION NUMBERS 
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Table 3.2 

Values Used in S/N Equations 



Table 3.1 Equation s d 

System # Designation ♦ Used Pe/sec HI N Pe/sec 
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Table 3.2 (Cont'd) 



The bandwidth terms were obtained using the following equations 


Af = Ax Ay A).Tf (3. 10) 

2 K 

where A^ = Number of pixels in x direction 

Ay = Number of pixels in y direction 
A>= Number of spectral regions 

Tf = Frame time (1/sec) 
n^= Number of detectors 
K = Scan efficiency 

Notice in Figure 3-1 that under low background conditions the 
systems are essentially signal noise limited and the point 
scanning systems have the better S/N. However, under high 
background conditions the multiplex systems are superior. 

3. 3. 4 Source Photoelectron Value Set Constant 

Figure 3-2 is a replot of Figure 3-1 where N^ is held constant 
at a nominal value and varies. All the other values in the 
equations are the same as in Table 3-2. From these curves it 
is obvious that the type of system to be selected is heavily 
weighted by the expected background term). 

4 Candidate System Trade Off Matrix 

A tradeoff matrix for the twelve candidate system is shown in 
Table 3. 3. An explanation of the table heading follows: 

a. S/N - This column shows the calculated S/N ratio 
for each system for a constant target flux (Ng) of 
5 X 1Q4 pe/sec and a low background level (Nj^) 
of 2 x 103 pe/sec on top and a high background 
level (Nq) of 1.95 x 10^ pe/sec below. 
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NUMBER OF PHOTOELECTRONS FROM SUM OF DARK CURRENT, ELECTRONS, X-RAYS & PROTONS 
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Figure 3-2. Signal- to- Noise Ratio vs Number of Photoelectrons From Additive Noise (Nq) For A 
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b, S/N improvement - setting system number 12 (the present 
BUV instrument) at 1, a calculation was performed 
to show the improvements achieved in the other system. 
The number in the top is for a low background and 
the number in the bottom is for a high background. 

c. Number of detectors - This column contains the 
number of detectors used in each system. 

d. Field of view - This parameter is the total field of 
view that the system images without scanning. 

e. Dwell Time/FOV - The time allowed for one field 
of view. 

f. Bandwidth - The system bandwidth. 

g. Space Radiation - This parameter is a calculated 
weighting factor derived from the S/N ratio under high 
background conditions. The factor is normalized to the 
Number 1 system. Increases in the value mean that 
the system is degraded by background. 

h. Optical System Field of View - This parameter implies 
the complexity of the objective lens. 

1. Image ^^otion Compensation - Does the system require 

image motion compensation? Since the vehicle is moving 
does the scene move during a frame. Y is yes and 
N is no. 

j. Sun Angle - Is the system effected by sun angle? Y 
is yes, N is no, 

k. Statistical Independence - Is each sample statistically 
independent"? ie. Does the MTF have a quasi square 
roll off or a 3db per octive roll off"? Y is yes, N is no. 
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stray Light - Can the stray light be controlled in the 
system - U is unknown, N is no. 

m. Number of Multiplexing Masks - The number of masks 
required and the % transmission of the masks is contained 
in this column. 

n. Comments - Pertinent comments on the system. 


From Table 3. 3 the following reduction in candidate systems 
have been performed: ' 

a. Systems 1 and 2 can be eliminated since system 1 
has a problem with the monochromator (see comments) 
and the dwell time is too long. System 2 is eliminated 
Decause the mask probably cannot be made. This 
system may also have monochromator problems as 

in system 1. 

b. System 4 is eliminated because the S/N improve- 
ment is not adequate. 

c. Systems 7, 8, 9, and 10 are eliminated because they 
are essentially the same as systems 3, 4, 5, and 6, 
besides which they require an object plane scanner. 

d. Systems 11 and 12 are eliminated because the S/N 
improvements are not adequate. 

Thus we are left with systems 3, 5 and 6. If the background is 

n 

high (2 x 10 pe/sec) then system 3 is a viable solution. How- 
ever, since the predicated mission does not have backgrounds in 
this region, system 3 was also eliminated. 
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Final Candidate Systems 


3. 5 

By process of elimination the selection of only two systems for 
further analysis was made. The two systems are systems 5 
and 6. Both of these systems perform the spatial scanning 
with a 15 element photomultiplier tube array in the direction 
perpendicular to the flight path and by the forward motion of 
the vehicle in the direction of flight. This type of system is 
called a pushbroom scanner. In the spectral domain the 
scanning is accomplished by either a Hadamard mask or in the 
same manner as the present BUV system. 
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PERFORMANCE EVALUATION 
General 


4. 0 
4. 1 

The two candidate systems selected in the candidate tradeoff 
analysis (Section 3. 0) are a pushbroom spatial scanner usirig 
either convention or Hadamard spectral scanning. For the 
conventional spectral scanning system, the performance has 
been based on an extrapolation of the BUV sensor, and, therefore, 
does not have to be proven. On the other hand, the Hadamard 
spectral scanning has inadequate performance history and thus 
theoretical and experimental results are presented to help 
characterize performance. Since the spatial scanning uses the 
pushbroom technique, which is well understood, only the 
Hadamard Multiplex Spectrometer will be discussed. 

4. 2 Hadamard Multiplex Spectrometer 

4. 2. 1 General Description 

A multiplex spectrometer is sensitive to many of the total 
number of individual wavelengths in the spectrum simultaneously. 
In comparison, a conventional scanning spectrometer (like the 
BUV) is only sensitive to a single wavelength at any one time. 
Thus, the multiplex spectrometer is more efficient. Mechanically, 
the single exit slit of a conventional spectrometer is replaced 
by a multislit mask which allows more light to pass through the 
spectrometer to the detector. 

Also, Hadamard spectrometers are transform spectrometers, 
since the raw data stream from the instrument is a mathematical 
transform of the optical input spectrum. Hence, they share with 
Fourier-transform spectrometers the requirement for postmeasure- 
ment decoding by a digital computer to recover the optical 
spectrum. The theory of operation of Hadamard-transform 
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spectrometers is covered in Appendix A. If one expresses the 
basic measurement operation in matrix notation, several things 
become apparent. First, that the multiplexing operation is 
essentially one of optically "writing a set of simultaneous 
linear equations in the n unknown spectral power densities" 

(that is, the unknown power of each of the n discretely resolved 
wavelengths in the spectrum). Next, that the "coefficient 
matrix" of this set of simultaneous equations simply describes 
whether a given wavelength is transmitted or blocked for any 
given measurement. Hence, the coefficient matrix is binary , 
corresponding to the two possible transmission states (say "1" 
for a transmitted wavelength and "o" for a blocked wavelength). 

The mathematical requirement for a solution to this set of 
equations to exist is that the coefficient matrix be "orthogonal", 
hence we are dealing with "binary orthogonal matrices". By 
definition, these are "Hadamard matrices." (the name' for any 
binary orthogonal matrix), or are derived from them, hence the name 
of the technique. 

The codes based on the Hadamard matrix can be demonstrated to 
be optimum when a single Hadamard mask is used for spectral 
multiplexing. (The Fourier transform, for instance, is theoret- 
ically a factor of two worse than the cyclic Hadamard-core 
"S-code" commonly used for HTS spectrometry (Tai and Harwit, 
1976 ). 

It should be pointed out that the Hadamard transform is a spatial 
transform, not a frequency or time-domain transform. Essentially, 
the Hadamard transform multiplexes a brightness distribution in 
space (the spectrum dispersed at the exit focal plane of the 
"disperser"), and assumes that the brightness distribution is 
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constant with time. As such, it is more analogus to a photogra- 
phic plate than to a time-resolved photomultiplier tube. 

4. 3 Sources of Noise: Effect on Transform Spectrometers 

Multiplex spectrometry should only be used if the resulting 
spectrum has an improved signal-to-noise ratio (SNR), or if 
greater spectral coverage or resolving power can be attained. If 
such advantages do not accrue, multiplex techniques generally 
imply increased expense and should be avoided. 

The discussion of SNR advantages is complex, because many 
different factors can contribute to the noise of a spectrum. 

In this section we will list some of the sources of noise and 
show how they scale. This should give a reasonably clearcut 
view of the conditions which make spectromodulation desirable 
and those which do not. 


4.3.1 Sources of Noise 

A. Photon Noise 

The radiation incident on a detector can be considered to 

consist of photons arriving randomly distributed in time. 

Strictly speaking, this complete randomness is true only of 

coherent radiation at low intersity levels, but for most 

purposes, that restriction is not important. If N^ photons 

interact with a detector in a given time interval, the random 

1/2 

signal fluctuations will be of order N^ , and the ratio of 

1/2 1/2 

the signal to photon noise is ~N^/N^ = N^^ , 


N ~ 
P 



(4. 1) 
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B. Background Noise 

If the detector sees a considerable background signal, a 
further contribution to the noise can result. If the back- 
ground - say thermal emission from the detector housing, or 
from the spectrometer - is unmodulated, it merely adds on to 
the photon noise due to the radiation whose spectrum is to be 
analyzed. If the radiation again provides photons in a 
time interval, and the background produces N2 photons 
which interact with the detector, the total noise will be 

V 

Nb~(N2+N^) 2, (4.2) 

The background noise however need not always be uniform. 

A distant source whose spectrum is analyzed through a 
changing atmosphere may produce detector readings which 
reflect time-varying atmospheric emission. 

The atmospheric noise N then has a form 

atm 

where a(t) is an effective variability in the atmospheric 
emissivity. Its absolute value ranges from 0 to 1. Its 
effects on the readings can strongly depend on the modulation 
employed in obtaining readings, since variations in 
atmospheric emissivity have peaks at certain modulation 
frequencies. Unfortunately, such variations appear to be 
wavelength dependent, and not enough is known about this 
type of noise to guide, say an astronomer wishing to make 
infrared observations at the most effective modulation 
frequency. 

4-4 
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C. Modulation or Scintillation Noise 

A source of noise similar to that given by Eq. (4. 3) is the 
unintended modulation the incident beam undergoes on its 
way to the detector. Sources of this noise may be atmos- 
pheric scintillation, or atmospheric patchiness - variable 
cloud cover, for example. Effectively the noise then has a 
character 


N ~/3(l)N,. (4. 4) 

md 1 

and again the various Fourier components of the time varying 
function )3(t) generally are not known. 


D. Detector Noise 

Many detectors have intrinsic noise levels much higher than 

the photon noise. This noise often is modulation frequency 

dependent and the experimenter normally elects to make 

his measurements in a frequency interval in which the noise 

is low. Once this has been done, the rms noise integrated 

1/2 

over an observing time T is proportional to T .If the 
noise in unit time interval is n^: 


^d 


V 

n^T^Z. 


(4. 5) 


For some types of detector, the noise is also proportional 
to the square root of the detector area A. For these 
detectors 


N 


V 

A~"2^ 



where is the noise per unit detector area. 


(4. 6) 
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other detectors effectively have high noise at high modulation 
frequency. This comes about because some very low noise 
detectors have very high impedance, and a detector's 
response time constant may become so long that neither the 
signal nor the noise is properly registered in available 
measuring time intervals. In such a situation, amplifier 
noise will become smaller at higher modulating frequencies. 

As the observing time decreases, the noise effectively goes 
up, if the signal is maintained constant by increasing the gain: 

N ~T ^2. (4. 7) 

amp 

I 

E. Minimum Detector Size 

Sometimes the minimum available detector size is much 

larger than the finest resolution element of available masks. 

To make optimum use of the sensitive detector area, encoding 

may then be indicated. For minimum available detector size, 

this noise N . scales like N , or N , Eqs. (4. 5) and (4. 7) 
size d amp 

If the detector is larger than the minimum size, and falls into 

the N., (below) class, N . will scale like N N 
A H size A, H 

to the point where the minimum detector size is reached, and 

like N,, N , below, 
d amp 

F. Housing Noise 

When a detector is kept in a housing which radiates at 

wavelengths to which the detector is sensitive, the detector 

may be more strongly limited by this source of noise, N , 

H 

than by any other. When possible, the housing should 
then be cooled, but in some applications that choice may 
not be available, and this type of photon noise will then 
dominate. 
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4 . 3 . 2 


Optimum Operation 


Let us now examine the best operating mode given the various 
sources of noise just mentioned. Amplifier noise does not guide 

us towards or away from spectral multiplexing. In general, a 
spectrum having m distinct elements requires m distinct measure- 
ments, whether we measure one spectral element at a time or 
all of them simultaneously. Either way therefore, a maximum 
time T may be available for any one measurement of the entire 
spectrum, and this will fix the noise level. It is clear, however, 
that in the presence of this type of noise, only one spectrum 
should be obtained in the available time interval. If two indepen- 
dent spectra of the source are taken in the available time, the 

signal-to-noise ratio of each one will be lowered by a factor 

1/2 

of 2, and the sum of the spectra will be worse by 2 than a 
single spectrum taken in the same time interval. This happens 
because the noise components of the summed spectra grow in 
random walk fashion while the signals add linearly. This argu- 
ment, therefore, speaks against rapid scan spectra when the 

limiting noise is from the amplifier N . In most other situations 

amp 

rapid scanning provides a succession of spectra which can help 
to average out effects of atmospheric and/or modulation 
noises. Since multiplexed spectrometers gather more light onto 
a detector in all cases, their use is still to be preferred when 
amplifier noise is present. 

Let us next turn to the modulation noise component due to 
scintillation. Much of this noise can be avoided if a "total 
intensity" reference detector is used. In fact, this should always 
be done if modulation noise is suspected. The ratios of the 
instantaneous intensity measured by the reference detector should 
then be used as the primary spectral data. The only time that 


4-7 



such a procedure would not work is if the modulation is also 
wavelength dependent in an unpredictable way. Spectral 
multiplexing and high throughput are definite advantages if the 
ratio technique is employed. Otherwise, however, no advantage 
is gained by multiplexing. It is difficult to compensate for the 
atmospheric emission noise unless all signal time variations 
are attributed to this noise source and the reference detector 
signal is primarily used in a subtractive mode. Clearly, 
ambiguities arise when both atmospheric and modulation noise 
sources are present. If atmospheric noise is the dominant 
source of uncompensated noise, multiplexing will provide no 
advantage, and an unnecessarily high throughput will actually 
provide a disadvantage, since the noise may grow even faster 


than the signal. 
If photon noise 


dominates, Fellgett's advantage cannot be 


realized. By looking at m spectral elements simultaneously, the 

1/2 

noise of each measurement increases as m . The total time T 

available for looking at each spectral element increases as m 

and since the noise components of successive measurements 

1/2 

add randomly, the growth of noise will be proportional to T 

The total noise, at the end of time T, therefore, is proportional 
1/2 1/2 

to m T = m. The signal strength, integrated over this 
time interval, is also increased by a factor of m, and there is no 
net gain in signal-to-noise ratio through multiplexing. 


Detector noise N- always makes multiplexing worthwhile. Let the 
o 

spectrum consist of m spectral elements and the noise for any one 
measurement lasting time T/m be n^. Let the signal received in 
the same time interval be s. If each spectral element is seen for 
a total number of times ^ m, the signal at the end of the total 
observing period T will integrate to Since the noise is 
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present even when no signal falls onto the detector, the noise n^^ 

must be summed in random walk fashion over all m time intervals 

1/2 1 1/2 
and integrates up to n^^m . The SNR therfore improves by 

Since many detectors show increasing noise with increasing size, 

we must still examine how to optimally operate when area noise 

dominates. In this case multiplexing must often be accompanied by 

2 

'dedispersion", Fellgett (1958) , Dedispersion consists in 

passing the radiation back through the spectrometer, once the 

9 

modulation at the exit slit has taken place. Decker (1971) has 
shown how this can be readily accomplished in practice. 
Dedispersion can gather all the radiation back onto an area as 
small as the entrance slit of the instrument and in principle no 

increase in detector size is then implied by multiplexing. Just 

as for detector noise, Fellgett's advantage holds in 

this situation. If the radiation is not dedispersed, a bigger 

1/2 

detector is needed. The noise increases as A , the effective 
time increase for viewing the signal goes up by ^ A and the 
SNR decreased by ^ ~ it actually is lowered! 


4.3 .3 


Computation and Construction Tolerances 


Two final characteristics which must be included in any compari - 
son of spectrometric instruments are the computational complexity 
required to produce a usable spectrum and the construction 
tolerances necessary to insure correct operation. Both of these 
requirements are at a minimum for conventional scanning-mono- 
chrometer spectrometers: there is no necessary post-measure- 

ment computation whatever, and, with the exception of the initial 
manufacture of the diffraction grating, the constructional tolerances 
are merely those normally associated with the construction of 
conventional optical systems. 
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In the area of computational complexity, the Hadamard-transform 
spectrometer has a clear advantage over other multiplex systems. 
In principle, the Fast Hadamard Transform runs approximately an 
order of magnitude faster than the equivalent Fast Fourier Trans- 
form, per transform dimension, Pratt et. al. , (1969). When one 
takes into account the fact that Hadamard spectrometers typically 
operate with much smaller multiplex numbers than do Fourier 
spectrometers, and hence require the decoding of a much smaller 
data matrix, Hadamard-transform spectrometers typically require 
far less in terms of computer power and computer time than do 
Fourier-transform spectrometers. For instance, a typical FTS 
system might involve a minicomputer with 16k to 32k of core 
memory and also a high-speed disc memory containing up to 
1 million additional words; a typical Hadamard system might 
involve either a stripped-down minicomputer or a microprocessor, 

either with only 4k to 7k of core memory (and no disc or tape 
memory at all). In addition, the Hadamard transform lends 
itself to inversion by simple, special-purpose, real-time "hard- 
wired" computers or systems based on microprocessor chips and 
read-only memories, Slingerland, (1975). 

Similarly, in the area of mechanical complexity and construction 
tolerances, the Hadamard spectrometer shows a clear advantage 
over other multiplex spectrometric systems. Fourier spectro- 
meters, as they are interferometric instruments, require not only 
fraction-of-wavelength-of-light accuracy in the motion of the 
variable-length arm of the interferometer, but for satisfactory 
operation, also require the use of up to 3 separate interferometers 
operating simultaneously: the signal (infrared) interferometer, 
a reference interferometer looking at a helium-neon laser (to 
calibrate the interferometer fringes and synchronize the data 
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system), and a "white-light" interferometer (to provide proper 
phasing of the "zero order fringe"). Hadamard spectrometers, by 
contrast, require only conventional optical tolerances, except 
for the diffraction grating and the multislit code mask - which, 
of course, need only be manufactured once. Typically, then, HTS 
systems require construction and operation tolerances of the 
order of fractions of millimeters, while FTS systems require 
tolerances of the order of fractions of micrometers. 

4. 3.4 Echo-Correction 

Closely related to the matter of construction tolerances is the 
matter of error susceptibility and error correction capability. 

All multiplex systems are in one respect improvements over con- 
ventional scanning spectrometers, in that partial data loss 
produces only degredation of the overall signal-to-noise ratio of 

the spectrum, while similar partial data loss with scanning 
spectrometers result in "gaps" in the spectrum - total loss of 
information about the particular spectral region under study 
during the time of the data loss. The other side of that coin, as 
it were, is that all multiplex systems distribute isolated high- 
noise "bursts" throughout the entire spectrum, degrading its 
signal-to-noise, while monochrometers inherently isolate noise 
bursts to the specific spectral regions during which they occurred, 
rendering those areas essentially illegible, but preserving the 
signal-to-noise of the remainder of the spectrum. 

One can further ask what types of operational or manufacturing 
errors are most probable for a given type of system, what their 
effects would be on the resulting spectra, and what (if anything) 
can be done about the situation. The most probable error to 
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which the Hadamard spectrometer is subject would be one in 
the construction or motion of the encoding mask. Small random 
errors will, to first order, merely degrade the final spectral 
signal-to-noise ratio, and are not susceptable to systematic 
correction. 

Systematic errors in the construction (or motion) of the multislit 
coding mask will, however, result in a systematic pattern of 
"echos" in the output spectrum, (Tai, et. al. , (1975b). Tai and 
his co-workers showed that for masks whose transparent slits 
are systematically wider or narrower than they should be, a 
single input spectral line will be distorted to an output spectrum 
containing that "correct" line plus a set of small "echos". When 
the transparent slits are too wide, these echos are of equal 
height, and generally consist of one pair surrounding the line and 
a second pair bunched together and displaced a definite specific 
distance from it; when the slits are too narrow, both sets of 
echos still have the same amplitude, but the displaced echos 
are negative. The location of the echos is uniquely determined 
by the code design of the mask, and hence is the same for all 
masks of a given design. The amplitude and sign of the echos; 
their position is determined by the mask code itself! For 
instance, in the common case of an HTS mask of 255 slits, an 
error in the slit width shows up as a pair of echos spaced 24 and 
25 spectral channels from the input line; if the slits are too wide, 
all four echos are "positive" (same as correct line), if too 
narrow, the echos displaced 24 and 25 channels are negative. 

As the effect of error in the coding mask can be uniquely related 
to the single-spectral-line response of the spectrometer, it is 
thus possible to correct for the effects of this error in the process 
of the computer decoding; Tai et. al, (1975b) have devised such 
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a correction procedure; This procedure consists of essentially 
three steps: First, a single-line spectrum is run with the code 
mask under consideration, to determine the amplitude and sign 
of any errors present, and the appropriate scaling constant of the 

echos is calculated; this need only be done once during the 

/ 

initial system calibration for any instrument. Next, the decoding 
program is modified so that it constructs a second "replica" of 
each (arbitrary) transformed spectrum, where the "replica" has 

been multiplied by the echo scaling factor determined during 
initial calibration. This scaled replica is then "cycled" through 
the spectrum and added (or subtracted, depending on the sign 
of the mask error) from the original spectrum so that each spectral 
point has had added (or subtracted) from it the scaled values 
corresponding to the known echo locations. This sounds compli- 
cated, but can, in fact, be accomplished with 4 or 5 additional 
lines in the FORTRAN "decoding" program. 

Figure 4-1 shows the effect of this correction on the spectrum of 
an 11. 018 m line-filter taken with a 255-slot Hadamard-transform 
spectrometer. The mask error was approximately 5% of the width 
of an individual slit, or about 0. 008 mm, and the resulting echos 
have been corrected to far less than the spectral noise level. 

This capability to correct for the dominant system error as part 
of the routine decoding procedure is unique to Hadamard 
spectrometry; Fourier spectrometry errors are seldom correct- 
able and typically result in depredation of spectral quality. 

4. 3. 5 Spectral Manipulation 

We would like to add a final note on one "multiplex advantage" 
which Hadamard shares with all other multiplex spectral tech- 
niques: the ease of post-measurement spectral manipulation. 
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Fig. 4-1. Mask-Error-Ecno correction of HTS spectra (11. 018|* line-filter source). 
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due to the necessary use of a digital computer in the measurement 
cycle. In both Fourier and Hadamard spectrometry, the output 
spectrum exists in the computer memory as an array of numbers, 
and hence is available for ensemble spectral averaging, 
normalization, ratioing, feature subtraction, mixture reduction 
or any other mathematical manipulation desired by the experi- 
menter. This has been described, among others, by 
Griffiths (1974) and Hirschfeld (1976) and is one of the major 
reasons for the growing acceptance of Fourier techniques among 
analytical spectrochemists. 

We wish to note here that these manipulative possibilities are 
fully applicable to Hadamard spectrometry. As one example, we 
refer to Figure 4-2: the upper trace (D4) is a single-beam spectrum 
of polystyrene film, seen in the region from 7. 5 to 15^, which is 
seen superimposed upon a source spectrum (R9, central trace) 
containing, among other things, a sharp order-sorting-filter 
"cut-on" and a rather deep "resonance" from this filter. The 
lower trace (D4/R9) shows the effect of computer ratioing of the 
sample spectra with the stored IR source "reference" spectrum 
to produce a "corrected" sample spectrum. 

4. 4 Hadamard Mask Design 

4. 4. 1 Mask Description 

A standard Hadamard mask utilizes a cyclic code which makes 
it possible to construct a single double-length multislit mask 
which contains all of the patterns necessary to encode the 
spectrum. Since, in the HROM instrument, the spectral bands 
of interest are not contiguous and are not equally spaced, a 
standard cyclic Hadamard mask cannot be used. Instead, a 
discrete field stop which contains the monochrometer exit 
slits that define the required twelve spectral bands, must be 
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Fig. 4-2 Compensction of sample spectra (here, Polystyrene film) by 
computer ratioing to stored reference spectra (here, IR 
glower soijrce). 



used in conjunction with a Hadamard mask. The field stop is 
shown in Figure 4-3. Following the field stop is the Hadamard 
mask. This mask moves normal to the dispersion direction 
and uses a separate pattern for each measurement. Thus, 
the total mask is composed of n individual patterns. The slits 
on the mask line up with the exit slits in the field stop as shown 
in Figure 4-3. Use of the field stop reduces the accuracy 
requirements on the Hadamard mask, since the spectral bands 
are defined by the field stop. 

4. 4. 2 Selected Code 

The code selected for this application is a 15 slot cyclic code, 
having the form: 

10001 0011010111 

The complete code for all 15 positions is shown in Figure 4-4. 
Since the code has 15 positions, and the instrument only 
requires 12, an inefficiency is present. To remove this in- 
efficiency, three imaginary spectral bands are measured 
^3’ ^14’ ^15' T^sse three bands equal 0. However, 

by using these bands, 15 measurements of 15 unknowns are 
performed, and the inefficiency is removed 

4. 4. 3 Mask Size 

The radius of the Hadamard mask is determined from the following 
formula: 


R 


d 


1. 2 H 


s 


2 tan 360° 
2n 


+ A 

s 

A7 H 

s 


(4.8) 
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10001 00110 10111 
00010 01101 01111 
00100 11010 11110 
01001 10101 11100 
10011 01011 11000 

00110 10111 10001 
01101 01111 00010 
11010 11110 00100 
10101 11100 01001 
01011 11000 10011 

10111 10001 00110 
01111 00010 01101 
11110 00100 11010 
11100 01001 10101 
11000 10011 01011 


Fig. 4. 4 


Hadamard Code 



where: 


H = The height of the slits 
s 

1. 2 = Factor which makes slits 20% oversized 

n = Number of mask positions = 16 
= Focal plane width 

s 

A y 


Ak 


Total spectral range = 3398 ^ - 2555 X = 
843 ^ 


AY = Spectral resolution =10^ 

W = Width of slits = A 



2 

A = Area of slit = 0. 153 cm 


s 

The proposed mask contains 15 positions for the Hadamard code 
and one blank position. The blank position is used to measure 
a dark noise every frame. Thus, the number of mask positions 
is 16. In section 4. 5 the area of the slits A required for a 

s 2 

30/1 signal-to-noise ratio is defined as 0. 153 cm . Therefore, 
substituting values into the equation, we get: 


R 


d 


3. 0164 H 

s 


+ 


12. 8979 
H 


s 


(4.9) 


This equation is tabulated in Table 4-1. Where is the slit 

heights, W^ is the slit widths, L is the width of the focal 

plane and R, is the radius of the Hadamard mask wheel. The 
d 

equation is also plotted on Figure 4-5. Notice that there is an 
optimum value for a given area slit. For the proposed mask, the 
slit height is 2. 068 cm, the slit width is 0. 074 cm and the 
radius of the wheel is 12.475 cm. 
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4. 5 Performance Comparison 

Since the conventional spectral scanner is based on the BUV 
instrument and the Hadamard spectral scanner has just been 
shown to be possible, it remains to compare the two sensors 
from a performance viewpoint. The two systems will be compared 
again using signal-to-noise ratio. However, this time the values 
used in the equations will be based on realistically derived 
system parameters. 


4. 5. 1 Number of Photoelectrons in the Signal (Ng) 

The number of photoelectrons in the signal can be calculated 
from the formula: 


Ng = NAg fi Tq A X T) 10‘ 


(4. 10) 


where N = Scene radiance - given in work statement as 

o 

2.5 X 10*’'^ erg/cm2 - sr - A - s 


Ag = Area of slits =2 .137 cm high by 0,072 cm wide 

^ = Solid angular field of view = 7° x 7° = 

0.01493 sr (given in work statement) 

Av o 

= Spectral resolution = lOA (given in work 

statement) 

Tq = Optical efficiency = 2% (same as present BUV) 

o 

X = Wavelength of interest = 2555A = 0.2555/xni 

o 

V = Quantum efficiency of PMT@ 2555A = 15% 
(PMT catalogue valve) 

n 

10 = Conversion factor ergs to watt. sec. 

h = Planck’s Constant = 6.625 x 10 sec^ 

c = Speed of light = 2 .99793 X 10^^ M m/sec 
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Substituting the values into equation (4. 10) we get: 

^ (2. 5x10"'^) (0.154) (0.0149) (0.02) (10) (0.2555H0. 15) (10~^) 
s (6.625 x10“34)(2.99793x10^'^) 

Ng = 2.21 x10^ photoelectrons/sec. 

4. 5. 2 Bandwidth (a£j^) 

The noise bandwidth for the two systems can be calculated from 
the following expression: 


where 


= Qx Qy e X Ft 

^ 2nd K u y (>> X 


(4. 11) 




w 


X 


Total field of view in the x direction = 104° 

Total field of view in the y direction = 7° 

o o 

Total spectral range 2555 to 3398A = 843A 

Instantaneous field of view in the x direc- 
tion = 7° 


W . 



W y 

e. 


w ■ 


K 


= Instantaneous field of view in the y direc- 
tion = 7° 

o 

= Spectral resolution = lOA for each of 12 bands 

= 15 resolution elements (pixels) per line 

= 1 resolution elemt high 

= 16 spectral bands (bands 3 ^ 13 , 3 ^ 14 / and 

a, 15 = o ) for the Hadamard system and 
12 spectral bands for the conventional system 

= Number of frames per second = V7.38 sec. 

= Number of photomultiplier channels = 15 

= Scan efficiency - 98% (0.01 sec per 22.5® 
step for Hadamard mask) 75% (littrow type 
mirror for conventional system). 
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substituting the values into equation 5.2 we get: 


Afn = as) (1) (16) 

H 2 (15) (.98) (7.38) 


Afnjj = 1.106 Hz For the Hadamard System 

= (15) (1) ( 12 ) 

2 (15) (.75) (7.38) 


Afn^, = 1.084 Hz For the Conventional System 

4. 5. 3 Signal -To -Noise Ratio fS/ N) 

A. Hadamard System (s/n)H 

The signal-to-noise -ratio for a Hadamard System was derived in 
section 3.3.2 and is: 



where 


Ns (MP) 
2N Afnjj 


1/2 


N 


1/2 


sMP 

2N 


nl/2 



(4. 12 ) 


(S/N)jj = 


N 


s 


P 


M 

N 

AfnH = 


N 


D 


Signal-to-noise ratio for a Hadamard 
System = 30 (given in work statement) 
Number of photoelectrons in the signal = 
2.21 X 10 pe/sec (calculated in Section 
5.1.1) 

Number of spectral bands = 16 (12 from 
work statement, X i 3 » x 14 and =0 ) 

Number of spatial resolution elements = 1 
(15 PMT's, one pixel per tube) 

Number of Hadamard masks = 1 
Noise bandwidth for a Hadamard System = 
1.106 Hz (calculated in Section 5.1.2) 
Number of photoelectrons from additive 
noise = ? 
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Solving equation (4. 12)for and substituting values of 1 for M 
and N. The equation becomes: 


Nd = NsP 


N. 


-1 


2tV 

WH 


Substituting values into equation (4. 13)we get: 


Nn = 2210(16) 1^ 2210 

2 2(30)^ (1.106) 


•l" 


Nd - 1.947 X 10'^ photoelectrons/second. 


(4.13) 


1 ) 


N«>>N, 


For noise in signal limited operation equation (4. 12) reduces to: 


/S\ ^ ^ 

\N /H 1^2 Afnpj 


1/2 

1/2 


(4. 14) 


2 Afnpj 

and substituting values into equation (4. 14) we get: 

(n)h 


2 ) 


N, < Ni 


For additive noise limited operation equation (4. 12)does not change. 

4 

Substituting values into the equation and letting Nj^ = 1.62 x 10 
pe/sec the signal-to -noise ratio becomes: 


(n)h 


= 22 .83 > N 
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3) Ng = Nd 

When the additive noise in signal are equal equation (4. 12) does 

not change. Substituting values into equation (4. 12) and letting 

= N the signal -to -noise ratio becomes: 

D s 



29. 80 N 

s 



B. Conventional System (S/n)c 

The slgnal-to-nolse ratio equation for a convention system is: 



N 
s 


^ 1/2 


(4. 15) 


where the symbols are the same as for the Hadamard System and 
the values are the same except for ^fn^ which equals 1.084 Hz. 


Solving equation (4. 15)for we get: 


= N 
D s 


N 




r-J ^ Afn 


-1 


(4. 16) 


substituting values into equation (4. 16) we get: 


= 2.796 X 10 pe/sec 


1) Ns>>Nj3 

For noise in signal limited operation equation (4. 15) reduces to: 

1/2 


= N 

iN/c s 


— 1/2 

Afn 


(4. 17) 


Susstituting values into equation (4. 17) we get: 


(n) 


= 45. 15 when Ng>>NQ 


2) Ns< Nd 

For the additive noise limited case-equation (4. 15) does not change. 
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4 

Substituting values into the equation and letting N^) = 1.62 x 10 
pe/sec the signal-to-noise ratio reduces to: 



15.64 For Nj^> Ng 


3) Ng=Nc 

When the additive noise and the noise in signal are equal, equa- 
tion (4. 15)does not change. Substituting values into equation 
(4. 15) and setting Ng = Nq the signal-to-noise ratio becomes: 

= 31.93 For Ng = Nq 
C. Crossover Point 

At what value is the signal-to-noise ratio for the Hadamard System 
equal to that of the convention i’ystems ? This value can be found 
by equating equations: (4. 12) and (4. 15) as follows: 
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rearranging equation (4. 21) we get: 
nl/2 - 


2 FA fnjj 

N3P/2 

1/2 

Afnc 1/2 j 

fWs +Nd 

L J 

1/2 


(4. 22) 


and 

r.1/2 


[ 




Afn^ , 2 FAfnH^^YNsP 


D 


1/2 


(4. 23) 


DJ c ti p 

squaring both sides of equation (4. 23) gives: 

pTn^ Afn^ = 4p2 Afnjj [^N^P/j + Np] (4.24) 

Performing the multiplication indicated in equation (4. 24) yield: 
PAfn^ Ng + P Afn^ Nq = 4p2Afnjj NgP ^ 4F2Afnjj (4* 25) 


Collecting like terms in equation (4. 25) we get: 
2F^Afnjj NgP - PAfn^ Ng = PAfn^ Nq - 4F^ Afn^j 
Factoring equation (4. 2^ gives: 


NsP 


Solving equation (4. 2 7) for we get: 


2F^Afnjj- Afn^ 


PAfn^ - 4F^Afnj^ 

- 


_ 


2 


Nd = N^P |2F^AfnH- Afn^ 




'PAfn^ - 4F^Afnjj j 


(4. 26) 


(4.27) 


(4. 28) 


If we let 

P = 15, Afnjj = 1.106 H^, A fn^ = 1.084 and Ng = 2210 pe/sec 
and substitute these values into equation (4. 28) we get: 

- 1.084 

(4. 29) 


35360 


2.212F" 


t 


17.344 - 4.424 F 


I ” 
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Finally since (|) = (|)^ ' f = 1- 

Np =3.09xl0%e/seo (|jjj = ^ 

4. 5. 4 Conclusion 

Figure 4-6 shows the signal-to-nolse ratio for a Hadamard 
System and a conventional system as a function of 
additive noise. Notice that the conventional system starts 
out with a better S/N ratio however when the additive noise 
level reaches 3090 pe/sec. , the two are equal and if 
more noise is added the Hadamard System excels. Again, 
as in Section 3, the signal-to-noise ratio cannot be used 
directly to select a final candidate system. 
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SELECTION OF RECOMMENDED SYSTEM 


5. 0 

5. 1 System Comparison 

Having completed the performance evaluation (Section 4. 0) and 
shown that either system can be characterized, it now remains to 
determine which system is preferred, both from considerations of 
performance and engineering simplicity. 

5. 1. 1 Performance 

Since the only difference in the two systems is the spectral 
scanning, most of the performance criteria can be met equally 
by both systems. These include such items as precision, 
dynamic range, calibration, linearity, fidelity and thermal 
influence. From a signal-to-noise ratio standpoint, either 
system can meet the 30 to 1 requirement. However, the 
conventional system is superior when noise-in-signal limited 
and the Hadamard is superior when additive noise limited. 
However, calculations show that the additive noise will be 
above the crossover point. Therefore, the Hadamard system 
is slightly superior. 

5. 1. 2 Engineering Simplicity 

Both systems use light pipes and 15 photomultipliers, 15 ampli- 
fiers, 15 A to D converters and 15 processing electronics; thus, 
electronically they are similar. From a telemetry standpoint, 
the Hadamard code has an advantage, since dropped bits or 
words only manifest themselves as a general lowering of the 
signal-to-noise ratio of the decoded scene, not as loss of 
information. Also, the simple rotary motion of the encoding 
disc is easier to mechanize than the oscillating motion required 
for conventional spectral scanning. 
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5. 1. 3 Recommended System 

The recommended system is the pushbroom scanner for spatial 
scanning and the Hadamard Scanner for spatial scanning. The 
selection was made on the basis of a slight performance 
advantage and a more substantial advantage in mechanization. 
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6. 0 SYSTEM MECHANIZATION 

6. 1 System Description 

The design of the HROM sensor is based on the experience 
gained from the BUV instrument flown on NIMBUS 4. An 
outline drawing of the HROM instrument is shown on figure 6. 1. 
Basically, the HROM system is composed of the entrance 
optics which define the instrument field of view and direct the 
incoming radiation onto the first collimating mirror. The 
collimated energy is directed to the first grating where the 
radiation is dispersed. The dispersed energy is then recollimated 
and dispersed again by a second grating. The doubly dispersed 
radiation is then encoded by the Hadamard spectral mask, 
collected by light pipes and transmitted to the photomultipliers. 
The light pipes collect the spectrally encoded radiation in such 
a manner that the spatial orientation is preserved. Outputs 
from the photomultipliers go to separate signal processing 
electronics and on to logarithmic A/D converters that interface 
with the spacecraft telemetry. 

6. 2 Optical Design 

The basic design of the optical system is very similar to that 
of the BUV instrument flown on NIMBUS 4 and the SBUV scheduled 
for NIMBUS G and is shown in Figure 6. 2. 

The BUV spectrometer is an F/5 double Ebert-Fastie mono- 
chromator using two 2400 1/mm gratings and two 250 mm focal 
length spherical collimating mirrors. Double dispersion is 
achieved at the exit slits by using three-corner-mirror transfer 
optics between the two monochromators, which makes the 
dispersion of the second monochromator additive. The use of 
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Figure 6. 1 HROM Instrument 
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double dispersion with a double monochromator system permits 
the entrance and exit slit widths to be twice as wide with twice 
the radiant flux throughput as for a single monochromator with 
the same spectral bandpass. In addition, the magnitude of 
stray light rejection is squared for the double monochromator. 

For example, if the stray light level for a single monochromator 

-3 

is 10 of the input radiant flux, then a double monochromator 

_ g 

of the same design should have a stray light level of 10 of 
the same input radiant flux. The entrance and exit slits are 
25 mm tall and 1. 5 mm wide and produce a 10 X nominal 
spectral bandwidth. 

The BUV monochromator uses no fore-optics thus its FOV is the 
11, 5° square angular acceptance of the F/5 monochromator. 

A calcite Lyot type depolarizer is used to minimize the sensitivity 
of the instrument to the polarization state of the backscattered 
ultraviolet radiation. The fundamental differences between the 
HROM spectrometer and the BUV spectrometer are these: 

1. The HROM spectrometer will use permanently fixed 
gratings and the full spectral image surface from 2500 % 
to 3400 ^ for the full slit height will be scanned by a 
Hadamard mask. 

2. An objective lens is required in front of the entrance slit 
to image a 104° FOV on the 20 mm entrance slit height 
and a 7° FOV on the 10 ^ equivalent entrance slit width. 

3. The astigmatic spectral image surface for defining image 
details along the entrance slit height is located about 

10 mm beyond the best astigmatic spectral image surface 
for the Hadamard mask use. This astigmatic entrance 
slit height imaging surface is divided into 15 equal zones 
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which add up to equal the full height of the entrance 
slit image. Each of these 15 zones had a width which 
encompasses all the spectral radiant flux passing through 
the full 2500 A to 3400 ^ spectral image width, 10 mm in 
front of it. The spectral radiant flux which passes 
through each of the 15 zones is then conducted to a 
separate photomultiplier cathode. 

6. 2. 1 Monochromator 

The HROM monochromator has an optical layout similar to the 
SBUV with mirror transfer optics. The three-mirror transfer 
system has the advantage of not causing monochromator 
defocusing as a function of wavelength. One mirror is mounted 
singly while the other two mirrors form a roof mirror system with 
two precision prisms mounted together. This type of block 
construction should enable the roof mirrors to be fabricated very 
accurately, and then to be mounted together accurately and 
securely. This block approach also minimizes the volume 
required for mounting. While the gratings for the SBUV Subsystem 
monochromator will be made from a different master and will 
probably have a different ghost pattern, the SBUV slit transfer 
function should be very similar to that of the BUV instrument 
for most wavelengths over the 1600 ^ to 4000 % wavelength 
range. The slit transfer function of the SBUV may be slightly 
different near 1600 R or 4000 A due to slit curvature effects. 

If 2400 1/mm fixed gratings are used the spectral image width 
is 2 inches at the intermediate slit and 4 inches at the exit slit. 
This means the monochromator entrance/ exit slit center line, 
and intermediate slit center line must be 4 inches further apart 
than on the BUV, and the 15 photomultipliers must receive spectral 
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radiant flux from a 4-inch spectral image width. 


1200 1/mm gratings reduce these spectral image widths by half. 

600 1/mm gratings quarter them, so the intermediate slit spectral 
image width is 1/2 inch while the exit slit spectral image is 
1 inch. The monochromator which has been layed out for this 
study uses 600 1/mm gratings which gives a nearly square 
spectral image with a spectral image width of roughly 2 0 mm and 
a height equivalent to the entrance slit height of 20 mm. The ray 
traces for the monochromator are shown in Figures 6. 3, 6. 4 and 
6. 5. 

The spacing between the entrance/ exit slit center line and the 
intermediate slit center line must be increased about 25 mm. 

The excellent stray light rejection capability of the BUV instru- 
ment implies that the HROM stray light rejection capability should 
be excellent, since the same stray light control techniques will 
be used. This should allow accurate photometric measurements 
of weak spectral regions in the presence of strong solar 
emission lines or continuum. The response of the BUV instrument 

looking at the zenith sky from ground level is shown in Figure 6. 6. 

o 

The stray light level at wavelengths below 2800 A is six decades 

o 

below the solar continuum at 3200 A. The intensity falloff which 

o o 

occurs between 32 00 A and 3400 A is due to the reflective 
attenuation by the coating on the collimating mirrors. 

The possibility that stray light due to multiple diffraction by the 
grating may pass out through the exit slit of the HROM mono- 
chromator was studied. The result of this study once again 
confirms that all wavelengths which are multiple-diffracted by 
the first grating either go out the entrance slit or impinge upon 
the monochromator walls on the entrance slit side of the 
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Figure 
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Fig. 6. 6 Response of the BUV Instrument looking at 
the Zenith Sky from ground level. 
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monochromator. 

6. 2. 2 Objective Lens 

The objective lens must have a short focal length to produce 
a 52° half angular FOV at the upper or lower ends of the 20 mm 
tall entrance slit. This lens is used at only F/5 but it must have 
a very large diameter to properly focus the monochromator 
acceptance rays for the full 20 mm tall entrance slit. 

Fused silica has too low an index of refraction to be used for 
the positive element of this lens, thus sapphire is used. 

For an inline objective lens system, the 104° FOV, unfortunately, is 
vignetted by the photomultiple mounting system, therefore, the 
entrance optics have been rotated 90° by a fused silica prism. The 
prism surfaces have been given negative correction to reduce 
chromatic aberrations. 

6. 2. 3 Detector Optics 

The assembly drawing Figure 6. 2, shows an arrangement of 15 
photomultiplier tubes relative to the spectral image surface. Fifteen 
identical fused silica plates are used to conduct the spectral radiant 
flux from the 15 zones of the spectral image to separate photomulti- 
pliers. 

These plates are each 1.33 mm thick, and plane parallel in a side 
view. Seen in a plan view, they are slightly trapazoidal tapering 
from a 20 mm width at the spectral image surface to an 8 mm width 
at the photomultiplier. 

6. 3 Detector Module 

The HROM has been designed to use the EMR Type 510 photomulti- 
pliers, which is physically much smaller than the photomultipliers 
previously used in the BUV and SBUV systems. The EMR Type 
510N-06-07 photomultiplier uses a bi-alkali photo- 
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cathode and a suprasil window. While these photomultipliers 
may be considered somewhat developmental they should have 
flight experience in time to be used for the HROM program. 

The HROM uses the type N bi-alkali photocathode which has a 

high (20% typical, not including window) and flat quantum 

o 

efficiency at wavelengths less than 4000 A. 

6. 3. 1 Window Selection 

The PMTs use fused silica as a window material. This material 
was chosen because it exhibits much less scintillation emission 
in the presence of Bremsstrahlung photons resulting from high- 
energy electrons. Experience with the NIMBUS-4 BUV, which 
employed two PMTs with sapphire windows, indicates that 
scintillation-produced background would interfere with the 
minimum signals to be processed by the HROM if sapphire 
windows were used. EMR PMTs are made in a modular manner, 
with dynode elements in the form of circular stamped plates 
separated by glass rings. PMTs are potted with dynode divider 
resistors in place. The dynode resistors are metal glaze types. 

6. 4 Mechanical Design 

In this design, all optical elements are stationary except for 
the diffuser assembly and the rotating Hadamard Mask. A 
layout drawing of the instrument is shown in Figure 6. 2. 

6. 4. 1 Mechanical Configuration 

The monochromator is constructed in a monolithic housing. Cast 
from A356 aluminum alloy, the housing is precision machined 
and provides direct mounting of all optical elements on lapped 
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integral seats. The detector and diffuser assemblies also 
mount directly on the housing. This direct mounting technique 
minimizes assembly time and assures position and alignment 
stability of all optical elements. To minimize the effects of 
thermal defocusing and thermal gradients within the telescope, 
the reflective elements are mono-metallic with the housing 
(all are aluminum). The differential thermal coefficient between 
the A356 alloy housing and the 6061 alloy mirrors is only 1 
microinch per inch per degree F. Using aluminum mirrors main- 
tains the geometrical relationships of the optics with the 
housings over a broad temperature range. This approach 
eliminates the need for an expensive structure to provide low 
thermal expansion and instead provides a high thermal 
diffusivity to prevent the occurrence of thermal gradients. 

Aluminum mirrors also reduce weight and simplify mirror mounting 
techniques. 

The system can be thermally controlled by passive techniques 
and thermally isolated from the Spacecraft mounting structure. 

The instrument utilizes a multi-layer aluminized mylar 
insulation blanket to minimize orbital thermal variations. An 
internal heater and a nadir-facing radiator may be included to 
accommodate various potential duty cycles and power dissipations. 
These techniques can maintain the instrument within a 10°C 
temperature range and provide instrument temperature stabilities 
in the order of 0. 5°C. 

The instrument is designed for a single-plane mount with the s/c. 
The weight of the proposed design is 15 kg and the volume 
is 0. 305 m high by 0. 305 m wide by 0. 61 m high. 
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The collimating mirrors and grating substrates are made of 
0. 2 inch thick, nickel plated, stabilized 6061-T6 aluminum. A 
rib pattern is added to the rear of the mirrors to facilitate 
fabrication to the required figure of one-fourth wave. 

The fundamental reasons for the use of the metal optics is to 
overcome tilt, de-center and defocusing errors as a result of 
differential thermal expansion between the optical parts and 
their housings. Metal optics are particularly suitable for 
optical systems which must withstand severe temperature 
changes and yet remain in alignment. 

6. 4. 2 Diffuser Plate Mechanism 

The diffuser plate is a transmissive element similar to the 

design employed on the BUV instrument flown on Explorer 55. 

The diffuser plate mechanism must deploy the diffuser plate into 
the field-of-view on command. It must also provide protected 
storage to prevent surface degradation due to long-term exposure. 

The diffuser plate is driven by a permanent magnet stepper- 
gearhead combination, and the entire mechanism is a removable 
assembly. 

Verification of diffuser plate position is accomplished with 
LED-photodiode pairs. Only two positions need be sensed — 
stowed and calibrate, so a simple pattern of small holes in the 
gear web will suffice. The rotating output member is balanced, 
so that motor detent torque multiplied by the gear ratio is 
sufficient to hold the output against rotation due to launch 
vibration. 
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6. 4. 3 


Hadamard Mask Drive Mechanism 


The Hadamard disc will be a 10 inch diameter thin (0. 001 inch) 
metal plate. It will be stepped in 16 equally spaced angular 
positions holding in each position approximately 400 milliseconds 
and stepping to the next position in 10 milliseconds. The step 
will be accomplished by means of a small constant speed DC 
motor dissipating 2 watts. The motor will drive a momentum 
storing wheel (a flywheel) which will have an escapement 
mechanism drawing power from the stored energy at the proper 
intervals. 

Since the energy required to accelerate/decelerate the Hadamard 
disc is only a fraction of a watt, the mechanical efficiency of 
the escapement need only be better than a few percent, which 
is easily obtainable. 

6. 5 Electronic System 

6. 5. 1 Signal Processing Electronics 

A block diagram of the signal processing electronics is shown 
in Figure 6. 7. The system employs 15 electrometer preamplifiers, 
one for each UV detector. The output of each detector is inte- 
grated for the mask dwell time (a little under 400 milliseconds) 
and converted to a 12 bit digital code. This is compatible with 
the required 1% radiometric precision as shown in Appendix B. 

If the output of the integrator at any code position exceeds 87. 5% 
of the full scale value of the ADC, a flag is set which is used to 
reduce the gain of the preamplifiers at the completion of the mask 
cycle. Conversely, if the output of the integrators at any mask 
position other than the opaque position fails to exceed 25% of 
the full scale value of the ADC, a flag is set which is used to 
increase the gain of the preamplifiers at the completion of the 
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Fig. 6-7 Signal Processing Electronics Block Diagram 








mask cycle. One mask cycle is completed in 7. 38 seconds, the 
time required for the scene viewed by the spectrometer to advance 
by one resolution element. During this mask cycle, 240 measure- 
ments are obtained, one for each detector for each of the 16 mask 
positions. 

One of the mask positions is completely opaque. It thus yields 
a black level calibration for the system and also provides 
verification of proper tracking between the mask cycle and the 
telemetry main frame timing. 

Because of the long duration of the measurement cycle, the 
possibility of commutating a single ADC among the 15 signal 
channels could be entertained. It seems, however, that with the 
advent of inexpensive microcircuit converters the resulting 
savings in hardware would be minimal, and would be more than 
offset by loss of the redundancy inherent in the use of indepen- 
dent parallel data channels. 

6. 5. 2 Hadamard Mask Driver 

The Hadamard Mask Driver advances the mask by one position 
after the integration of each sample. The time interval while the 
mask is in motion is not available for signal integration, and 
should therefore be minimized. The required cycle time for each 
mask position is a little less than 400 milliseconds. This is 
based on the spatial and spectral resolution requirements of the 
instrument. An advance time of 10 milliseconds would correspond 
to an equivalent collection efficiency of 98% and would therefore 
seem acceptable. 

To achieve such rapid motion with an acceptable power level, a 
low inertia system is required. The present primary candidate 
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system is a thin disk driven by a DC motor with a flywheel. 

Aside from the Hadamark Code slits, the disk contains a code 
pattern suitable for detecting and identifying sixteen positions 
by means of an electro-optical sensor. 

In normal operation, the disk is driven from each position until the 
next position is detected. One of the sixteen positions, 
corresponding to a fully opaque mask, is identified as the home 
position. At this position, a tracking counter is reset which is 
advanced each time the disk is advanced. The state of the 
counter is subsequently compared to the position indicator and 
sets an alarm flag if disk position and tracking counter indication 
ever fail to agree. This alarm flag is transmitted to the ground by 
way of telemetry to alert the data reduction team that the corres- 
ponding measurements are suspect. 

A Back-Up Operating mode, used in case of failure of the position 
sensing system, advances the code disk by a predetermined 
number of stepper motor steps to the next position. Ground- 
commanded single step capability allows proper phasing of the 
wheel in this mode. 

6. 5. 3 Calibration 

Spectral and intensity calibration cycles are provided. 

Spectral calibration is accomplished by a suitable spectral lamp 
illuminating the entrance aperture of the instrument. The elec- 
tronic system provides power to start and operate the spectral 
lamp at the required intensity. 

Intensity calibration is accomplished by inserting a sunlit 
diffuser in front of the mapper's entrance slit. The electronic 
system provides the electrical driving signal to actuate a suitable 
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mechanism. Electrical timing signals control power for the 
spectral lamp and for the diffuser mechanism in such a manner 
that a calibration cycle starts and terminates at the opaque 
position of the Hadamard mask. The calibration cycles are 
thereby constrained to last an integral number of full cycles. 

This minimizes lost data and facilitates data reduction. Calibra- 
tion cycles are identified by a suitable flag in the housekeeping 
data bit stream. 

Black level calibration is obtained once each Hadamard cycle, 
by letting one position of the mask disk be fully opaque. 

5 5^ 4 Housekeeping and Diagnostics 

The HROM provides housekeeping data to assist data reduction. 

This data will include commutated binary digital signals as well 

as uncommutated analog signals. Digital signals include: 

An indication that the mask wheel is tracking. 

An indication that the amplitude cal lamp is on. 

An indication that the amplitude cal diffuser is in. 

An indication of the gain state of each detector channel. 

An indication of the operating modes of the instrument 
selected by command. 

Analog signals include: 

Indications of the following four temperatures: 

Detector Assembly 

Code Disk Drive Motor 

Electronics Assembly 

Low Voltage Power Supply 

An indication of the several output voltages of the low 
voltage power supply. 

An indication of the output voltage of each high voltage 
power supply. 
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Digital housekeeping signals are obtained by sensors directly 
monitoring the condition of interest, or by sensing the power 
actually supplied to the actuator or exciter. 

Temperatures are measured by miniature thermistors fastened to 
the assembly or part to be monitored by means of a thermally 
conductive bonding material. The conductance of the thermistors 
is measured by conditioning circuits with a calibrated transfer 
characteristic. Suitable buffering is provided to prevent loading. 

Voltages are monitored by divider networks. When negative 
voltages are monitored, an operation amplifier provides inversion 
for telemetry compatibility. 

6. 5. 5 Commands 

The HROM contains provisions for receiving serial command data 
from the spacecraft command system. The instrument command 
circuitry comprises a command receiver, decoders, and control 
registers. The receiver accepts command data from the space- 
craft interface in serial form and provides a parallel output that 
is valid immediately after command data is received. The 
decoders route the appropriate command data to the proper control 
register and the control registers accept and store the command 
data and use it to control the operational deployment of the 
instrument. Control registers are generally updated at the 
end of a Hadamard code cycle so that instrument configuration 
is changed at the opaque position of the mask. 

Commands are used, among other functions, for the following 
purposes: 

To activate or deactivate the spectral calibration 
source 

To activate or deactivate the intensity calibration 
diffuser 
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To switch the mask disk drive from the manual mode 
to the back-up mode or vice-versa 


To enable or disable any one or more high voltage 
power supplies 

To advance the code disk by a single motor step in 
the back-up mode of operation 

If compatible with the spacecraft command system, commands 
will be encoded in such a manner that each of these functions 
can be controlled independently of the others without reference 
to the prior command configuration of the spectrometer. 
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APPENDIX A 


Theory of Hadamard Spectrometers 

A - 1 Codes For Hadamard Transform Spectrometry 

Now consider a Hadamard-transform spectrometer which is 
encoded only in the exit focal plane. 

Let W be an NxM matrix whose rows represent the M encoding 
masks successively placed in the exit plane of such a 
spectrometer. A very desirable characteristic of W is that it 

e + 

be cyclic; that is, the (i+1) row of W is obtained by shifting 
th 

the i’^" row cylically one place to the right. 

Such a mask has two advantages. First, it can be self- 
supporting and therefore permits the construction of a spectro- 
meter which requires no transmission materials. In operation, 
the mask is stepped one slit width along the length of the mask- 
i, e, , in the direction of dispersion - for each successive 
encoding position. 

More important, perhaps, cyclic codes enable us to avoid the 

tedious construction of N masks each having N slots, for a total 
2 

of N slots. Instead we construct only one mask with 2N-1 
slots. The cost of mask construction is reduced by N and 
the design of the advance mechanism is considerably simplified 
since the bulk weight of the masks also decreases as ~'|’N. 


J 
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- 2 Singly-Encoded Instruments 

Sloane, et. al. (1969) have developed the following theory: 


The optical spectrum whose shape is to be estimated is 

spatially dispersed and the band of interest partitioned into 

th 

N channels. The average energy in the j channel, after a 

selected integration time, is denoted bySf'. The measurement 

process consists of observing the spectrum through M masks, 

th N 

the energy through the i mask being where W, 

J = 1 ij j 1 

til 

(« 1 , . . . , is the i vector of mask elements. The 

photodetector adds a random noise d. to the signals;^ w, 

1 "-j = l ij J 

and yields a measurement 


N 

Zojij'l'j, 1 = 1, 2 M. (A-1) 

The noise u. has the following properties 0; is 
independent of the signal: <J^; successive measurement 
noises are assumed to be uncorrelated (<( 1 ^^ 0 if i / j) and 

<C ^represents an ensemble average. 

In order to extimate |'K| by we need at least as many 

measurements Mas there are unknowns N. Furthermore, at 

least N distinct masks are needed to estimate the apectral 

shape. Hence, assume vectors of observations ^ = (’T^i • • • , 
channel energies^ = . . . , '1'^^), measurement noises ii. = 

T T 

(•^f, . . . , and a matrix of masks W = (W^, . . . , W^) = 

('jj.). ^The T stands for transpose). With this notation ='j^W +£. 

A 

An estimate of^of^is a function of the observations, (’j), 
expected to lie close to As a measure of the accuracy of the 
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estimate we adopt the mean square error criterion: We minimize 
« = < ( 0 -^) 

For purposes of computational convenience we agree to restrict 

A A 

^ to be a linear function of the observations, = tjA, for some 
matrix A; in the absence of more detailed statistical knowledge 
concerning the anticipated spectral shape or the photodetector 
noise characteristics, this seems appropriate. Before an 
essentially uniquely best experimental design can be derived, 
however, we still need to make the assumption that the estimator 
is unbiased; i. e. ,<( This assumption can be defended on 

the following grounds: 1) we need an estimator which, on the 
average for a large number of applications, yields the true value; 
2) the unbiased estimator can be shown to be desirable when 
there is a large uncertainty in ^relative to our prior knowledge of 
the spectrum and the measurement noise power. 

We must now select those matrices A, W that minimize c subject 
to the constraint we have set. First we see that 

<^^>= <^A>= <»7>A. (A-2) 

Furthermore, 

<V>= <'^>WK«'>= <^>W, (a-3) 

Hence, unbiasedness requires that WA = I, the identity matrix. 

If M = N, W is square, and A = W \ If W is not square but is 
N X M (M > N), we can use the generalized inverse. 

A = W^iWW^) 

Having solved for A in terms of W we now select W to minimize e 
<(M) (^-^)^> . We note that 
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(A-5) 


\{aWA+ ‘'A-^=»' W^(W w"^) ^ 

By assumption <^t>^vy = ^i. Thus 

e/cr^= Trace [w'^(WW^) '\ww^) "V]. (A-6) 

2 -1 -1 T 

If M = N, this simplifies toe/a = Trace (W (W ) ). The 
optimum experimental design will be completed if we can find 

2 

that W - subject toa>.. = 0 or 1 or -1- that has the minimum c/o- . 

ij 

A - 3 Choice of Masks 

We just showed that in the case M = N, when there are as many 
measurements as unknowns, the matrix W = (w. .) of mask elements 

ij 

-1 -LT 

should be chosen so that = 0 or 1 or -1, Tr (W (W ) ) 

is as small as possible. Three possible choices for the matrix 
W are given here. See Appendix A of Sloane et. al . (1969) for 
the characteristics of Hadamard matrices needed to understand 
the following results: 

If H is an N x M normalized Hadamard matrix, and G is in the 
(N - 1) X (M - 1) matrix obtained by deleting the first row and 
column of H, and S is obtained from G by replacing + l's by O's 
and - I's by + I's then three possible choices for the matrix W 

of mask weights, together with the corresponding values of 
Tr (W ^(W ^)^), are as follows: (Table A - 1 ) 

We may note that 

T 

(a) Another choice for W is R , where R is the matrix 

obtained from G by replacing - I's by O's. This gives 

a trace, however, which is slightly but uniformly worse 
T 

than that from S . 

(b) All the matrices given in Table A-1 are superior to the 
single slit case (where W = I). In the single slit case 
Trace = N. 
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TABLE A - 1 

Mask Weights 


Matrix W 

Trace 

Elements 


1 

+ 1. -1 


2-2/N 

+ 1. -1 

s'^ 

4-8/N+4/N^ 

+ 1, 0 
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(c) Nelson and Fredman (1970) , Philips (1972) , and 

Sloane and Harwit (1976) have shown that the Hadamard 
matrix and its variations describe a set of optimum masks. 

A - 4 Optical Arrangement 

T 

A multislit spectrometer using the S code works as follows: 
Radiation which passes through the entrance aperture is 
rendered parallel and directed towards the dispersive element. 

The dispersed radiation is collimated and focuses upon the multi- 
slit mask at the exit plane of the instrument. 

The radiation transmitted by the mask passes through the post- 
optics of the system and impinges on the detector. We then 
obtain a spectrum by sequentially stepping M masks at the exit 
plane and recording the detector output for each mask. The 
inversion procedure described above then enables us to recover 
the spectrum. 

In some instances, depending upon instrumental design and the 
choice of code, it may be advantageous to make use of both the 
reflected and transmitted radiation of a mask. This scheme 

T 

would be particularly useful for masks which utilize the H and 
T 

G matrices. Here the + I’s would represent reflecting slots 

T 

and the - I's would represent transmitting slots. For the H 

matrix all elements of the first column would be + I's and the 

T 

masks corresponding to each row of H would reflect the first 
spectral element at all times. The remaining N -1 slots would 
be stepped in the usual fashion. Such instruments have not yet 
been constructed. 

Infrared spectral measurements normally require that the radiation 
be chopped. One can then realize codes with O's and I's by 
chopping between the transmitted radiation and a standard 
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source. The codes with + I's and - I's could be realized by 
chopping between transmitted and reflected radiation. 

A - 5 Comparisons with Other Dispersive Spectrometers 

To make a fair comparison, we assume: (1) a fixed total 
measuring time T; (2) constant slit widths; and (3) an equal 
number N of unknown spectral elements to be estimated. For the 
doubly encoded system there are initially 2N - 1 unknowns 
!)’■■■ ’^0’ * ■ ■ ’'^N - r comparison 

we will suppose we only wish to estimate the N central elements 
^(N-l)/2 ^o. . ^ (N-D/2 


(taking N odd for convenience). Of course, we can still obtain 
some estimate of the ends, but these errors will not appear 
in the comparisons. 


As a measure of performance we take, as before, the total mean 
square error in all the unknowns: 


2 

^ total 


(N+D/2 

E 


t= -(N+D/2 



(A - 7) 


Table a- 2 compares three different grating spectrometers. The 

first column is for a single entrance and exit slot. N measure- 

2 

ments are made in time T, with a mean square error o’ in each, 
and the second column is for a singly multiplexed instrument 
with an exit mask. 
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TABLE A-2 


Comparison of total mean square error for two grating spectro- 
meters in estimating N unknowns 


N 

Conventional 

Singly Encoded 

3 

3.2 

2. 25a^ 

7 

1 ? 

3. 06d^ 

11 

\\? 

3. 36d^ 

19 

19d^ 

3. 61a^ 

N 

n/ 

|2 - 2/(N + 1)1 ^0^ 


(ff is the mean square noise in a single measurement made in 
2 

time T/N; o' = (constant) N/T). 






APPENDIX B 

Dynamic Range Requirements for High Resolution Ozone Mapper 


Statement of the Problem : 

It is desired to calculate the dynamic range required to transmit the 
Hadamard coded signal of the High Resolution Ozone Mapper without 
significant degradation due to the noise (or digitization level) of the trans- 
mission channel. The High Resolution Ozone Mapper is a 16 - channel 
spectrometer. Twelve spectral data channels are coded by Hadamard 
techniques . 


Discussion: 

The proposed instrument achieves the 12 elements of spectral information 
by 16 code mask positions. Hence the analysis for a 16 element instrument 
is appropriate. We therefore assume 


n = 16 spectral elements 

n = 16 mask positions (i.e., 1 Hadamard sample = 1 signal 
measurement per mask position) 
n = 16 time elements 

N Photons incident, total spectrum, total scan 


During each time element, N/n photons are incident on the instrument and, 
since n/2 mask elements are opaque, on the average the detector sees 
N/2n photons, N/n^ from each transparent mask element. 

The incremental signal for one spectral element, i.e., the amount by 
which the signal from that element deviates from the average, is computed by 
subtracting the sum of the signals for the n/2 mask positions for which that 
spectral element is opaque from the sum of the signals for the n/2 mask 
positions for which the same spectral element is exposed. 


If this computed signal expressed as a number of photons detected) 

is generated by an actual differential in photons ( AN) associated with that 
spectral element, then A<J) and AN are related by 


A = ■^ mask positions 


^ mask positions 


/N an \ 
\2n n / 


— photons 
2n 


photons - 


AN 

2 


photons 


The noise or uncertainty in the measurement is then obtained as the 
quadrature sum of the noise associated with each time element. The noise 
for each of the n/2 elements for which the spectral element is exposed is 


Ee = 




— 4 . AN 
2n n 
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and the noise for each of the 
opaque is ^ 



n/2 elements for which the spectral element is 


Hence the total noise is 
E 




( N +AN) 


For the spectral element of interest, the total signal is given by adding the 
incremental signal A<l>= AN/2 to the average signal (per spectral element ) 

N/2n. 


The noise in this total signal is given by the quadrature sum of the noise 
components. The noise component for the incremental signa l is given by Erj, 
abov e. T he noise component for the average signal is'iN/ 2i of which 
l/nVN/2 is associated with each spectral element. The signal-to-noise 
ratio for each spectral element is then given by 

MN) 

and if n »1 and N»AN this reduces to: 


fl /i ^ 

V2 In N i 


SNRj^ = 

If we now add the effects of the noise (£) associated with the transmission 
of each Hadamard sample, the total noise (E^') becomes 


F ' — 


= (N +AN + 2nE^) 


and a noise l/n'^N/2 + nS^ is associated with the average signal for each 
spectral element. The total signal-to-noise ratio for each spectral element 


snr; = 


ji^N+AN 


and if n»l, N>>AN 



SNR^' = 1 


NI ^ 


and hence 




= n [l 

f N/n + A N N 
1 2 SNR/ ) 

1 - f] 
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but the highest signal (S) for one interval is 


and if N » A N 
S = 


and 


then 


£2 = 


N 

. AN 

2n 

T 

n 

N 


2n 


1 

/ 2 s + An 

n 

V 2 SNRx’ 


/ ^ \ 

1 

( 1 + 2 S I 

n 

V SNRx' / 


)■- i 


therefore if n = 16, SNR > 100, <-i2S 


S " 100 N 16 




1 


400 


or, in words, a transmission channel signal-to-noise ratio of 400 -r 1 is 
required for a radiometric accuracy of 1%. It can also be seen that the 
required signal-to-noise ratio is reduced when A N is not small 
compared to 2S, i.e., when the energy is not equally distributed among 
spectral elements. It can therefore be concluded that the required 400 1 

signal-to-noise ratio refers to a worst case situation. 


It should be noted that this signal-to-noise ratio is required if one 
measures the transform signal with respect to the dark level, and 
therefore the range of the encoder must extend from the dark level 
to the signal level. In several applications the required dynamic range does 
not extend down to the dark level. In an instrument that uses Hadamard 
techniques to code image data, for example, all samples remain fairly close 
to the average. Unfortunately, that is not so in a spectrometer, since the 
full range of signal levels is required to encode a spectrum that consists 
entirely of one strong emission line. 
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